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Abstract 

Massive runaway positrons are generated by runaway electrons in tokamaks. The fate of these 
positrons encodes valuable information about the runaway dynamics. The phase space dynamics 
of a runaway position is investigated using a Lagrangian that incorporates the tokamak geometry, 
loop voltage, radiation and collisional effects. It is found numerically that runaway positrons will 
drift out of the plasma to annihilate on the first wall, with an in-plasma annihilation possibility 
less than 0.1%. The dynamics of runaway positrons provides signatures that can be observed as 
diagnostic tools. 
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Positron, the anti-particle of electron, is a rare species in the part of universe where we 
reside. Since last century, man-made positrons have been generated in accelerators for sci¬ 
entific research, in nuclear reactors as the byproducts, and applied in different fields, such as 
medicine and material detection l|-l3|. Recent researches indicate that tokamak, a magnetic 
confinement fusion energy device, may be the largest artificial positron factory in the world 
4, l5|. In large tokamaks like JET and JT-60U, above 10^^ positrons are generated in a 
post-disruption plasma by runaway electrons 6-[^|. The dynamics of these positrons after 
birth in tokamaks is a noteworthy question that may yield valuable information about the 
runaway dynamics and disruption process in tokamaks. What is the fate of these positrons? 
Will they annihilate inside the plasma or on the first wall of the vacuum chamber? Be¬ 
cause the annihilation probability dependents on the path and velocity, the fate of runaway 
positrons is determined by their phase space trajectories, which are strongly affected by 
the loop electric field, the helical magnetic field, and the collisional and radiation effects. 
Incorporating all these factors, we study the phase space dynamics of runaway positrons in 
tokamaks and predict their final fate. 

When relativistic runaway electrons interact with the thermal electrons and ions, 
positrons are produced in the pair production process. In tokamaks, because the energy 
of runaway electrons is typically above 10 MeV, most of their “children” positrons are born 
relativistic and can be accelerated to runaway velocity by the loop voltage. In the toroidal 
direction, the runaway positrons are accelerated along the toroidal electric field and slowed 
down by radiation loss and collisions. As the energy increase, strong synchrotron radiation 
and bremsstrahlung radiation begin to dominate in the drag force, which finally balances 
the loop electric field force. The collisional drag from the background plasma becomes small 
after runaway positrons gain high velocities. Nevertheless, the collisional effect offers a mo¬ 
mentum transfer mechanism between the parallel and the perpendicular momentum through 
the pitch-angle scattering. On the other hand, the projection of the positron gyrocenter tra¬ 
jectory onto the poloidal plane is not exactly located on a flux surface due to the geometric 
effect of the tokamak magnetic field. Numerical results reveal that the circular orbits of 
runway positrons in the poloidal plane drift toward or against the major radius direction 
Bji (see Fig{T]), under the influence of the loop electric field along the toroidal direction ^ 
and the helical background magnetic field. Due to the drift effect, the runaway positrons hit 
the first wall of the tokamak within about one hundred milliseconds with energy as large as 
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Figure 1: The circular concentric magnetic flux surfaces and the coordinate systems. 

150 MeV. After the phase space trajectories are known, the annihilation probabilities along 
the trajectories can be calculated. For a typical positron, it is found that the probability of 
annihilation inside the plasma is only about 0.1%. Essentially all positrons generated in the 
tokamak will annihilate on the first wall of the vacuum chamber. This fact suggests that the 
annihilation spectrum from the wall can be analyzed to infer the dynamics and distribution 
of the runaway positions. Admittedly, there are many other loss mechanisms for the posi¬ 


tions. For example, stochastic field lines 
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24| induced by MHD modes and ripple field 


25| can result in transport or orbit loss for energetic positrons. The orbit loss considered 


in the present study is a neoclassical effect due to the toroidal geometry, and is probably 
the slowest loss mechanism among all possible mechanisms. The in-plasma annihilation 
probability of 0.1% obtained in our calculation is thus its upper bound. 

We now present in detail the study of phase space dynamics of runaway positrons. For 
theoretical analysis and numerical simulation, the gyrocenter of runaway positrons is de¬ 
scribed by a Lagrangian, which incorporates the tokamak geometry, loop voltage, radiation 
and collisional effects. Parallel momentum, perpendicular momentum, annihilation rate, 
and the drift orbit in the poloidal plane are numerically calculated as functions of time. The 
potential of runaway positrons as a diagnostic tool is discussed at the end. 

As the anti-particle, the positron has the same rest mass as the electron, denoted as me, 
but opposite electric charge, denoted as e. So the dynamics of runaway positrons are similar 
to that of runaway electrons. The Lorentz factor for a positron with momentum p is 
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where c is the light speed in vacuum, is the momentum component parallel to the back¬ 
ground magnetic field, and p± is the perpendicular momentum. Its synchrotron radiation 
drag force takes the form 




sin ^0 


( 2 ) 


where Ug = e^/dTreorngC^ is the classical positron radius, Vg = pj^/eB is the positron gyro- 
radius, 9 is the pitch angle defined by sin 6 ' = p±/p, and Rq is the major radius of the 
tokamak. The bremsstrahlung drag force is 


Fb = (^g// + 1) 27 - ^ , (3) 

where Ug is the number density of the background plasma, Zf>ff is the effective ion charge 
factor. The collisional friction force is 


Fr.= 
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where In A is the Coulomb logarithm, which varies slowly with the plasma parameters. To 
theoretically investigate the dynamics of runaway positrons in tokamaks, we describe the 
dynamics of positrons by the following Lagrangian 26| . 


L = (eAo + eA/ -h eAefp\ + p\\h) ■ x - 


(5) 


Here, Aq is the vector potential of the background magnetic field satisfying B = V x Aq, 
A; is the vector potential of the loop electric field satisfying 
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A effw is the parallel component of the effective vector potential corresponding to the drag 
force. 
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and b is the unit vector along the magnetic field. The magnitude of the effective vector 
potential in the perpendicular direction is 
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Without loss of generality, we consider a tokamak magnetic field with circular concentric 
flux surfaces, 


R ' 

Accordingly we choose its vector potential to be 

Aq = 


2R 


2Rq 


The loop voltage is set to be 


R 


Bor ^ 
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In the right-handed coordinate system {R, z) (see FigH]), the Lagrangian takes the form 

L = pRR + Pii + PzZ-^rrieC^ , (12) 


where 
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Because of the toroidal symmetry, i.e. dL/d^ = 0, the effective toroidal momentum is 
conserved, i.e.. 


dL 

Pc = —^ = const. 


This invariance determines the evolution of p\\ as 
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Meanwhile, there exits another conserved quantity, the effective magnetic moment p, defined 
by 


P = 


{p± + Aeff±f 

2meB 


( 18 ) 


which determines the evolution of p_L. If neglecting the higher-order terms caused by the 
toroidal effect and the poloidal field, Eqs. flTT)) and flT8|) give 
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Figure 2: The evolution of the parallel momentum (a), the perpendicular momentum (b), and 
the annihilation rate (c) of runaway positrons in tokamaks with different initial values of {p\\,p±) 
(normalized by niec). The loop electric field is Ei = 3Vjm. 
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which are consistent with the momentum evolution equations in [12|, l27H29|. After substi¬ 
tuting Eqs. llTTll and (lT8l) into Eq. llT^ and dropping the term p^^, the toroidal symmetry 
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Figure 3: Snapshots, taken for every 0.03s, of the circular orbit of a runaway positron in the poloidal 
plane, with initial position R — Rq = 0.1m, z = Om and initial momenta py = BnieC and pj_ = nieC. 
The loop electric field is Ei = bVjm. The position of the first wall is indicated by the dashed 
vertical line. 

leads to the reduced Lagrangian in the {R, z) 2D space 

L = pr{R,z)R + P;,{R,z)z-H{R,z). (19) 

This is the procedure of Routh reduction. Then the runaway positron dynamics in the 2D 
conhguration space, i.e., the projection poloidal plane, is given by the Euler-Lagrangian 
equation, 


R = 


z = 


dH/dz 
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( 20 ) 


( 21 ) 


dpz/dR — dpR/dz 

Equations flTTl) . (ITSl) . (1201) and (1^ determines the dynamics of runway positrons in phase 
space. Given the dynamics in the momentum space, the in-plasma annihilation probability 
of a runaway positron can be calculated according to 


Ran = / neCran^(dr , 


(22) 


where the annihilation cross-section for the positron-electron reaction is 


7 











27r 


(23) 


vrr: 


^ an. 


1 + 7 




137p(l 


g-27r/137p^ • 


7" — i ' ■ -^72 — 1 

To guarantee the long term numerical accuracy and fidelity, we have adopted a varia¬ 
tional symplectic integrator 0, 1^ . which discretizes the Lagrangian directly, to carry out 
the numerical simulation. For conventional integration algorithms, such the Runge-Kutta 
method, numerical errors from different time-steps accumulate coherently and the error 
grows without bound for long-term simulations. In comparison, the variational symplectic 
integrator can globally bound the numerical error for all time-steps, and thus are superior 
for simulating the runaway dynamics, which often involves many hundreds of thousands of 
turns in the poloidal plane. 

For the present study, we use the parameters from EAST as a specific example 0|. We set 
Bq = 3.5T, Rq = 1.7m, Ug = and In A = 10. The loop electrical field is Ei = 3V/m, 

and the minor radius is a = 0.4m. In typical EAST operations, the safety factor q varies in 
the range of 1 to 3, while the profile of q is rather fiat in the core region. Thus we choose q = 2 
in the calculation. For runaway positrons with different initial parallel and perpendicular 
momentum, their momentum evolution are plotted in FiglS After 0.3s, all the parallel 
momenta reach a steady value around 300meC, as a result of the balance between the loop 
electric field acceleration and the drag force resistance. The perpendicular momenta also 
evolve towards a steady value, though a little slower, due to the balance between the radiation 
loss and the momentum transfer from parallel direction through the pitch angle scattering. 
The green curves show the dynamics of what might be called a "backward runaway" positron 
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34| . whose initial parallel momentum is opposite to the loop electric field. It undergoes 


a deceleration under in the toroidal direction at first. However, its parallel momentum 
reverses sign after 0.01s and then increases like forward runaway positrons with positive 
parallel momenta. Its perpendicular momentum has a drop in the deceleration phase and 
passes a point of inflexion at 0.01s. It can be seen that runaway positrons of widely different 
initial momentum nonetheless approach attractor curves in momentum space. The existence 
of attractor curves is made possible because of the dissipation introduced by radiation effects. 

The possibility of annihilation for runaway positrons with different initial momenta is 
depicted in Fig|2](c). Because the annihilation cross-section becomes very large when the 
positron moves slowly relative to the background electrons, the annihilation rate of the 
















positron with negative parallel momentum grows rapidly during its turning around pro¬ 
cess. Nevertheless, it escaped from the doom of annihilation for a non-relativistic positron. 
Though its parallel momentum goes through zero, its perpendicular momentum is still rel¬ 
ativistic, which ensures that it is not annihilated in the slowing-down process and can be 
accelerated to runaway in the backward direction. Overall, the annihilation probability is 
very low for runaway positrons. The probability is only about 0.1% after 0.4s, which is long 
enough for runaway positrons to escape the tokamak plasmas through the drift motion or 
other transport/loss mechanisms in the poloidal plane. 

The motion of a runaway positron in the poloidal plane with loop electric field Ei = bV/m 
is depicted in Fig 131 The loop electric field and the toroidal magnetic field are in the 
direction, and the poloidal magnetic field is in the Eq direction. The counterclockwise 
circular orbits are snapshotted for every 0.03s. The outward drift of the circular orbit to 
the Er direction is evident. The drift velocity slows down with the radiation resistance 
in the toroidal direction increases. The drift velocity can be roughly estimated as Vdr = 
q{Ei + Eeff\\)/BQ j3|- The runaway positron finally hits the wall before t = 0.15s, with 
an in-plasma annihilation probability less than 0.1%. This example illustrates the fact that 
most of the runaway positrons in tokamaks can drift out of the plasma to hit the wall before 
annihilation within the plasma. If the safety factor q = q{R, z) is a function of the spatial 
coordinates like in real tokamaks, the drift velocity will explicitly depend on the spatial 
location. However, the space-dependent safety factor also leads to similar runaway positron 
dynamics and the same main conclusions as the constant-g assumption. For situations with 
larger loop electric field, the runaway positrons will drift faster. If changing the helical 
direction of the magnetic field or the direction of the loop electric field, the positron will 
drift inwards, i.e., towards the negative er direction, to hit the inner wall. It is observed 
that the radius of circular orbit varies as the positron drifts along the e/j direction, especially 


in the later stage. As mentioned previously, the neoclassical orbi 


loss mechanisms for positrons. Since other loss mechanisms 
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loss is just one of many 


37| are in general 
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35 


faster, the in-plasma annihilation probability observed in experiments should be even less. 
Finally, because the orbit drift does not depend on the electric charge, note that runaway 
electrons, which will have their own signature (like visible damage to the wall), will strike 
exactly on the same side of the vacuum vessel. 

Early in 1986, Surko et ah proposed to diagnose the transport process by injecting 
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positrons into tokaniaks 
was also studied 


38|. The annihilation spectrum of positrons in thermal plasmas 


39| . Now it is clear that the large amounts of positrons produced by toka- 


maks themselves can be used as a diagnostic tool as well. Since most of the positrons are 
annihilated outside the plasma, the positron diagnostic in tokamaks cannot detect plasma 


properties directly as in PET (Positron Emission Tomography) jl|. However, positron diag¬ 
nostic is still possible using our knowledge of the dynamics of the runaway positrons. The 
annihilation spectrum, which can be recorded by the gamma spectrometer, is very charac¬ 
teristic and easily identihable. The intensity, breadth, and shift of the annihilation peak 
in the spectrum reflect the properties of the plasma that determines the runaway dynam¬ 
ics. Moreover, the time history of the positron annihilation is also an important indication 
to some events in tokamaks, such as a disruption or rf heating, which produce a burst of 
positron runaways at a specihc time. The annihilation locations on the hrst wall and the 
emission directions of the gamma ray provide information about the phase space coordi¬ 
nates of the runaway positrons at the end of their journey, from which we can also infer the 
trajectories of runaway electrons in the phase space. Theoretical and experimental studies 
on these topics will be reported in future publications. 
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